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Abstract. To describe a quantum system whose potential is divergent at one 
point, one must provide proper connection conditions for the wave functions at the 
singularity. Generalizing the scheme used for point interactions in one dimension, 
we present a set of connection conditions which are well-defined even if the wave 
functions and/or their derivatives are divergent at the singularity. Our generalized 
scheme covers the entire U{2) family of quantizations (self-adjoint Hamiltonians) 
admitted for the singular system. We use this scheme to examine the spectra of the 
Coulomb potential V{x) = — e^/|x| and the harmonic oscillator with square inverse 
potential V{x) = {muj'^/2)x^ + g/x^, and thereby provide a general perspective 
for these models which have previously been treated with restrictive connection 
conditions resulting in conflicting spectra. We further show that, for any parity 
invariant singular potentials V{—x) = V{x), the spectrum is determined solely by 
the eigenvalues of the characteristic matrix U G [/(2). 
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1. Introduction 



Quantum singularity is a source of interesting physics and, at the same time, confusion. 
Even in its simplest form as a point interaction — now realized approximately as quantum 
dots (see, e.g., [1]) — it provides unexpectedly rich quantum phenomena such as duality 
and anholonomy [2]. When it arises as a divergent point of an infinite potential wall, it 
may admit quantum tunnelling through the infinite wall and, in some cases, can lead to 
an exotic quantum caustic [3] . However, it also poses the problem in its own treatment in 
quantum mechanics. In fact, if we look back the history of the one dimensional Coulomb 
potential, V{x) = — e^/|a;|, for instance, we find persistent disagreement over the possible 
spectrum for nearly a half century [4] . A similar confusion can be found for a system with 
the square inverse potential, V{x) = g/x^ [5]. 

These confusing circumstances arise due to the ambiguity in choosing boundary (or 
connection) conditions at the singularity, which in mathematical terms corresponds to the 
choice of self-adjoint domains for the Hamiltonian operator, 

H ^ - — ^ + V(x). (1) 

It has been known that a point interaction on a one dimensional line admits a U (2) family 
of self-adjoint extensions for the Hamiltonian, and that these are characterized by distinct 
connection conditions [6]. For a point interaction occurring at a; = 0, the connection 
conditions can be given by 

(t/-/)* + zLo(C/ + /)*' = 0, (2) 

where U G ^7(2) is a matrix characterizing the self-adjoint extension, / is the identity 
matrix, and Lq 7^ is a constant with dimension of length [7,8,9]. ^' and are boundary 
vectors 

U(-o);' v-v^'M);' ^""^ 

defined from the boundary values V'(='=0) = hnix^io V'(^) of the wave function if; and 
its derivative = The problem, however, is that this prescription of connection 
conditions may not be directly applicable to singular potentials V{x), because then the 
boundary values i/^i^O) and/or i(;'{±0) may diverge at the singularity and, accordingly, 
the vectors and in (3) become ill-defined. 

One of the purposes of the present paper is to provide a scheme in which the connection 
conditions (2) become well-defined even for singular potentials, and thereby furnish a 
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general framework for studying singular systems on a line including those mentioned above. 
Typically, at the singularity such systems allow for two independent square integrable 
solutions X — for the eigenvalue equation Hip = Eip for any E, and (at most) only one 
square integrable solution at a; ^ ±00. Systems of this type are said to be in the limit-circle 
case at a; = and in the limit-point case at infinity [6]. At a limit-point singularity, no 
boundary condition is needed for ensuring the self-adjointness of the Hamiltonian [7] , while 
at a limit-circle singularity, some boundary conditions are necessary in order to specify a 
self-adjoint Hamiltonian from the family of possible self-adjoint domains. This family is 
U (2) for these systems, which follows, for example, from the fact that each of the negative 
and the positive half lines has one square integrable eigenmode for any nonreal eigenvalue 
E [6,7]. An essential step toward the generalization of the connection conditions for such 
systems consists of replacing the boundary values of the wave functions with corresponding 
Wronskians. This idea has been proposed [10] for systems on a half line with a singular 
endpoint for which a family of boundary conditions is assigned. Here, we extend this 
to the full line, where now the family is given by U{2), in such a way that the connection 
conditions (2) remain to be valid with modified boundary vectors. (For a different scheme 
of providing the domains of possible self-adjoint Hamiltonians, see [11].) 

For illustration, we employ our scheme to analyze two models, one with the Coulomb 
potential and the other with the harmonic oscillator with square inverse potential V{x) ~ 
(ma;^/2) + gjx^ . We shall see that the various different quantizations discussed previ- 
ously for those models arise at different choices of the matrix t/, and that the spectra are 
dependent on the choice of t/; in fact, this dependence has caused the confusion concern- 
ing the spectrum in the literature. Interestingly, the spectrum depends on two parameters 
(the eigenvalues of the matrix C/), not all the four oiU e C^(2), and this two-parameter 
dependence of the spectra is shared by any parity invariant potentials Y{—x) = V{x) with 
singularity of the kind just mentioned. More precisely, we find that for those systems the 
space of spectra is given by the Mobius strip U{1) x U{l)/Z2. 

The plan of this paper follows the line of arguments stated here, that is, we give the 
generalized connection conditions in Sect. 2 and thereby analyze the two models in Sect. 
3. The two-parameter dependence of the spectra is then established in Sect. 4, and finally 
Sect. 5 is devoted to summary and discussions. 
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2. Connection conditions 



To begin with, we provide tiie general connection conditions by extending the con- 
struction proposed in Ref. [10] (see also [12,13,14]) from the half line to the full line. Let the 
potential V possess a singularity at x = on the one dimensional line X.^ The potential is 
assumed to be in the limit-circle case at a; = from both sides and in the limit-point case 
at a; — > ±oo, and regular otherwise. We first consider the maximum domain c L^{^) 
on which the Hamiltonian H can be defined as a differential operator, 



(4) 



= {ip ^ L^(X) I ip and ip' are absolutely continuous on 

every finite subinterval of X\{0} , Hip e L'^{X) } . 

The Hamiltonian is not symmetric on J-', since, for (p, ip E J-', we have 

J^dx [<p*Hip - mr^p] = ^ mr,i^]+o - w[<p%^u) , (5) 

where W[(p* ,'ip]±o are the limiting values for a; ^ ±0 of the Wronskian 

W[(P*,iP]{x) = (P*{x)iP'{x)-(P*'{x)iP{x). (6) 

Here we have utilized the facts that VF[(/)*, i/j] vanishes for x ±00, since the infinites are 
limit-point [7], and that it is finite in the limits x ±0 even if the two functions (f){x), 
ip{x) are divergent. This latter can be shown as follows. For e > 0, we introduce the space 
of functions 

J^^ = {ip E L'^{0,e) \ ip and ip' are absolutely 

continuous on (0, e) , Hip e L^(0, e) } . 

Note that JF^ contains T as well as a wide range of other interesting functions as well, 
including all the cigenfunctions of the differential operator H, which are square integrable 
in any finite neighbourhood of the limit-circle singularity a; = but not necessarily on the 
whole line X. Now, for (p^ip e J^^ and < 5 < e, 

fdx [^H^ - (Hcpr^p] = ^ mr^^u - w[r,m ■ (s) 

Both terms of the r.h.s. are finite. The l.h.s. is also finite, even if we let 5 — > 0. Con- 
sequently, \im.s^Q W[(p* ^ip]5 = W[(p* ,'ip]+o is finite. As one can see, this property holds 
actually not only in but even in J^^- The finiteness of W[(p* ,'ip]-o is proved similarly, 
with the aid of the analogously introduced J^-e- 



^ We use the symbol X to stress that it is dimensionful in contrast to the dimensionless real line '. 
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Since the r.h.s. of (5) is generally nonvanishing, any self-adjoint domain T) for H must 
be such a subset of JF that the r.h.s. of (5) is zero for all functions within T). Now we show 
how to characterize the possible self-adjoint domains via a connection condition at x = 
of the form (2), where the boundary vectors ^' are appropriately generalized with the 
help of a basic set of energy eigenmodes. Let (/?*^*\ for i = 1, 2, be two independent, real 
eigenmodes with eigenvalue -E, 



(9) 



for X ^ 0. The actual value of E is unimportant for our purposes. We note that these 
eigenmodes may not be square integrable on the whole line and hence may not belong to 
but they necessarily belong to .Fe and T-e- Consequently, the complex column vectors 



(10) 



are well-defined for ip e since they are constructed from finite quantities. Further, 
observing that we can rewrite the Wronskian (6) as 



W[cf>*,iP] 



W[(t)\ ^(^^] W[^P, ^(2)] - W[(t)\ ^^2)] w[^p, ^(^)] 



(11) 



we can express the r.h.s. of (5) in terms of the boundary vectors for and ip simply as 



2m 



(12) 



where $ and are introduced for (j) analogously to (10) for ip. If 0, ip are in a self- 
adjoint domain T> then (12) must vanish. In particular, for (p = ip, this condition reads 
$t$' = $'t$^ which, under the notations 

*(^) = *±iLo*' (13) 

with an arbitrary nonzero constant Lq, is equivalent to the equality of the norms ||^'*^"'"-' || = 
ll^'^"-*!!. This shows that and ^f^~^ are in a relationship 



U e U{2) 



(14) 
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which is nothing but (2). Different states (f), ip E V have to share the same U so as to 
make (12) identically vanish: 

42mI/o 



^[$(+)^^(+)_(C/$(+))t([/^(+))] ^0. 



By an argument analogous to the case of the half line [10], it is not hard to show that 
the connection condition (2) restricts the space ^ to a domain V = Vu on which the 
Hamiltonian is not only symmetric but indeed self-adjoint. Since all different Us specify 
different self-adjoint domains Vu-, here we can see again that, as in the case of point 
interactions, the family of self-adjoint Hamiltonians H = Hu allowed on the line with 
potential V{x) possessing a limit-circle singularity and limit-point behavior for x ±cxo 
is given by U{2). U will be called the characteristic matrix oi the self-adjoint Hamiltonian 
Hu. 

Various subfamilies of U (2) can be defined analogously to the case of point interac- 
tions. For instance, the 'separated subfamily' Q,^ where no probability flow through a; = 
is allowed is characterized by those U which are diagonal. Indeed, for diagonal U the 
probability current 

^'^"^^ ^ 2^^^*^' " ^"^^ = 2^ (^[^*' V'^'^] -W[r, ^^'^]W[i^, ^^^^]) (16) 

is seen to vanish at x = 0, and diagonal Us provide the cases when the boundary condition 
(14) does not mix the +0 boundary values with the —0 ones. Hence we have O3 ~ 
U{1) xU (1) C ^7(2), which are, in other words, the cases where the system consists of two 
independent half line systems. 

If the domain of a self-adjoint Hamiltonian contains only functions that are regular at 
the singularity, then one may choose as reference modes any basis of independent solutions 
satisfying 

(^W(±0) = 0, (^W(±0) = 1, v^^\±^) = -l, v^^y{±<d) = Q. (17) 

Under this choice, we find that the boundary vectors (10) reduce to the conventional 
form (3), which shows that our connection conditions are a natural generalization of the 
conventional conditions. Once generalized, however, we recognize that the normalizations 
(17) are not at all essential in presenting the connection conditions (2) at the singularity to 
ensure the self-adjointness of the Hamiltonian. This in turn suggests that the characteristic 
matrix U characterizes the singularity only with respect to the reference modes chosen, 
and this fact has been implicit in the previous treatment for non-singular cases based on 
the normalizations (17) . 
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3. Two models with singular potential 



We now employ the scheme just presented to analyze the two models mentioned in 
the Introduction. 

(i) One dimensional hydrogen atom 

The first model is the one dimensional hydrogen atom, which is governed by the Coulomb 
potential, 

V(x) = (18) 

fI 

This system has a long history of research, dating back to Loudon [15] who first gave a set 
of bound state solutions {En < 0) for the Schrodinger equation, 

Hi^nix) = EnMx), (19) 

in terms of the Whittaker functions. The spectrum obtained in [15] is 

4 

TTie 

^" = -:;iTT' n=l,2,..., (20) 

where each level is doubly degenerate. The system has later been examined by a number of 
other groups to obtain different spectra due to different choices of the connection condition 
at the singularity x = (see, i.e., [16, 4] and references therein). The connection condition 
adopted originally in [15] is the Dirichlet condition '(/'(±0) = which corresponds to the 
Friedrichs extension of the Hamiltonian [17], but other extensions are equally possible as 
we shall now describe. 

To apply our scheme of connection conditions, we first recall that in terms of the 
variables 

^/-2mEn / -m 

the Schrodinger equation (19) becomes 

d^V'n /a 1 



n^-iJ^"""- '''' 

This is a special case of Whittaker's differential equation, whose two independent solutions 
are the regular Whittaker function, 

M.{z)=ze-^^F{l-a,2;z), (23) 
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where F[a, 7; z) is the confluent hypergeometric function, and the irregular one, 

(^) = {^^(1 - 2; ^) [in^ + ^^(1 - «) - V'(i) - ^(2)] 

«+Z_.H(r + l)! '"^ /• 

Here r(a;) is the Gamma function, ip{x) = ^lnr(a;) is the di-Gamma function, and 
r-i r ^ ^ 

-4, = E 



(24) 



n=0 



1 11 



n + 1 — a n + 1 n + 2 



From the asymptotic behavior of the two solutions, one finds that W^^i{z) is square 
integrable whereas M^^i{z) is not. Thus the bound state must be of the form, 

M^) = w^a.i (1^1) {Nr e{x) + Nl e{-x)} , (26) 

where ©(a;) is the Heaviside step function, and Nr and Nl are constants to be determined 
by the connection condition at a; = 0. Note that, since asymptotically 

W^a.i (^) = S^~ + z[\nz + ^{l-a)- ^{1) - + Oiz'lnz), (27) 

as 2; — > 0, the bound state ipn{x) has finite limits at a; ^ ±0 whereas the derivative V'n(^) 
diverges there. 

To see which bound states are actually allowed by the connection condition (2), let us 
first fix the reference modes conformity with (9). We choose them as 



1 

'{x) = - 

^^'\x)^-T{l-/3)Wf,^{2K\x\) 



= -_M^,. (2«|a;|) [6(0;) - Q{-x)] , ^^^^ 



with 

■\/—2mE / —m 



which are analogs of (21) with replaced by some arbitrary E < 0. With these, the 
boundary vectors (10) become finite as they ought to be, and they are proportional to 
each other, 

^ = a(^A, *' = e*, (30) 



where 

1 ^ 2me^ 



ln^-^(l-a)+^(l-/3) 



(31) 



r(l-a)' " 
With these, the connection conditions (2) read 

[{U-I) + iLo^{U + I)]^ = 0. (32) 

For i/jn to be a nontrivial bound state, we need 

det [U-I + iLoiiU + /)] = det [L* - / + iLoi^D + /)] = 0, (33) 

where is a diagonal matrix appearing in the standard decomposition, 

U = V-^DV, V e SU{2). (34) 

In terms of the parameterization, 

D={ J-l' ^±e[0,27r), (35) 



and 



we find that (33) is satisfied if 



Locot(^). (36) 



e = -J- or -J-. (37) 

Thus, given the singularity specified by U, we can determine the spectrum of the bound 
states as solutions of (37). We observe that the spectrum depends only on the two angles 
(6'_|_, 0-) in the diagonal part of U, that is, the two eigenvalues of U. Later we show that 
this is in fact the case for all parity invariant potentials V{—x) = V{x) sharing the same 
singular property considered here. 

A particularly simple spectrum is obtained at the angles (6'+,6'_) = (tt, tt), i.e., at 
U = —I. The connection condition (32) then implies \E' = 0, and therefore we need a — 0. 
From (31) we learn that this is fulfilled for ct = 1, 2, . . ., reproducing the spectrum (20) 
with double degeneracy (since and Nl are chosen freely modulo the normalization). 
Note that, since (f'^'^\±0) = and (f'^^^ (=tO) = 1, the condition \1/ = is equivalent to 
demanding tlj{±0) = 0, which is obtained by the so called Friedrichs extension discussed 
in [17]. We, however, emphasize that the Friedrichs extension is a special self-adjoint 
extension belonging to the separated subfamily Qs where the two half lines are physically 
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decoupled and none of the scale parameters L± in (36) appears in the spectrum, but any 
other extensions whose spectrum varies with the parameters {9j^, 9-) in U through L± are 
equally possible.^ 

We next turn to the scattering phenomena of the Coulomb potential under our general 
connection conditions. (These considerations will be valid for the repulsive Coulomb force 
as well, with setting < 0.) For this, in place of (29) we consider a positive scattering 
energy Ej^ > and use 

Then the two independent solutions are still given by the Whittaker functions with the 
new a = and z = 2ik\x\. One may choose the following real combinations for a set of 
two independent solutions: 

M^) = ^M,^^^ i2ik\x\) [Qix] - Qi-x)] , 

"^^J^ (39) 
M^) = - 2 {r(l - il) W,^,i{^ik\x\) + r(l + i7) W_,^^. {-2ik\x\)] . 

In passing, we note that the phase part of r(l + z7), i.e., r]o = argr(l + z7) is called 's-wave 
Coulomb phase shift' while its modulus can be evaluated as 



|r(i ± Z7)| = Vr(i + z7)r(i-z7) = y (40) 

The set (39) is chosen so that the solutions satisfy W[(f)'^^^ {x), (j)^'^\x)] = 1 and normalized 



as 



0i(±O) = O, 0;(±O) = 1, 02(±O) = -1, (41) 



whereas 02 (±0) are divergent and cannot be normalized. From (23) and (27) the asymp- 
totic forms of the solutions for \x\ oo are found to be 

et^ 1 

^^^^"^^ |r(l + i7) I fc ^^^^1 -^ln2k\x\ + rjo) [e{x) - e{-x)] , ^^^^ 
02 (a;) ~ ~|r(l + i^)\e~^^ cos {k\x\ — 7 ln2A;|a;| + 770)- 

The general solution for scattering states is a linear combination of the two solutions, 

il;{x) = {n^^Ui{x) + N^rU2{x)} Q{x) + {iVi'Vi(x) + N^^U2{x)}Qi-x). (43) 



2 

The four-parameter family of extensions for the one dimensional Coulomb system has been argued 
in [16] in a slightly different scheme. 
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Now, from the asymptotic behaviors (42) one deduces that the incoming wave from the 
left is the one that behaves for x ± cxd as 



V'(a;)~Te*('=^-^'°s2'=^) (a; ^ +00) 

^ ^i{kx+'yln{-2kx)) _^ ^^-i{kx+'yln{-2kx)) _^ — oo). 



(44) 



This corresponds to the choice 



AtW = z/ce-^^-''°|r(l + i7)|T, N^^^ = -^.^^—-^T, (45) 



and 



N^^^^ = -ike-^^-''^°\r{l + i-f)\{R-e^''''), = - .^'^ ''^\. {R + e^^^°). (46) 

|r(l + z7)| 

To implement the connection condition at a; = 0, we need to find a set of reference 
modes, (p^^^x) and (p^'^\x) satisfying (9). For this we shall use the same set (42) of the 
solutions, (p^'^\x) = (f)i{x) for z = 1, 2, with reference energy E and, accordingly, with the 
corresponding parameter ^{E) obtained from (38). The boundary vectors (10) can then 
be evaluated as 



with 

p = ^ {f{E) - f{Ek)} , m = {21n7(«) - V(l - ilis)) - + il{s))} • (48) 

h 

Plugging the vectors (47) into the connection condition (2), and solving for the scattering 
matrix, one obtains 

l)- ^U-I)-Z(U + I) l(U-I)+^'L4U + I)l{l), (49) 
where we have used 

oj = ke-^'^\T{l + i^)\'^ + ipe^^. (50) 

Note that the Friedrichs extension U = —I allows no transmission T = 0. In fact, as seen 
easily in (49), this is the case for any diagonal U, which is expected from the fact that 
those U belong to the separated subfamily ^3. 
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To render the scattering data (49) more explicit, we use the decomposition (34) with 
D given in (35) and V parametrized as 



y = e*2'^2e^^''% /ie[0,7r], e [0, 27r) . (51) 
Then, in terms of the scale parameters (36) and 

X± = arg(l + za;L±), (52) 

the scattering formula reads 

(T\ ^_^^{2^o+x++x-) ( isin(x+ -X-)sin/xe-^^ \ 
\RJ \cos{x+-X-)-ism{x+-X-)cosi^J ' 

We note that this outcome depends on the choice of the reference modes, not only on 
the matrix U. This is a consequence of the fact that the combination of U and the reference 
modes, not each, is essential for the determination of the singularity, as mentioned earlier. 
The ambiguity in the choice of the reference modes is described by the group SL{2, R) on 
account of the reality condition of the modes and the normalization in the Wronskian (9). 
Precisely which combinations of the parameters are physically important is an interesting 
question and will be discussed elsewhere. 

(ii) Harmonic oscillator with inverse square potential 

Our second model is the harmonic oscillator with inverse square potential, 

V{x) = —-x^+g^. (54) 

In contrast to the previous example, we consider the repulsive case ^ > to examine the 
positively divergent potential, and add the quadratic term to render the entire spectrum 
discrete. To comply with the condition that the singularity be of the limit-circle case at 
a; = 0, we confine ourselves to 

0<9<^, (85) 

for which the Hamiltonian admits a U (2) family of extensions. To solve the Schrodinger 
equation (19), let us set 
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for a; > 0, and use 



1 / 8mg 2 f^^-. 

so that (19) becomes 



^^(^) + (a + 1 - ^(^) - ^ (a + 1 - An) /„(^) =0, A„ = ^. (58) 

This is just the confluent hypergeometric differential equation, and hence the two inde- 
pendent solutions for (19) are 



ci = 1 + a, 



C2 = 1 - a. 



(59) 



Since the solution for a; < can be found by setting x — > —x in (59), the general solution 
for the bound state is given by 

M^) = [N^U^:H\x\) + N^U^^\\x\)Mx) + [Ni'U^:\\^\) + nPcI>^:\\x\)M-x), (60) 

(s) (s) 

where the constants Nj^ and A^^ ^^^^ be restricted by the connection condition. Note that 
(55) guarantees that both of the two solutions are square integrable near the singularity. 
The entire square integrability is then ensured if the solutions vanish sufficiently fast at the 
infinity a; — > ±00. From the asymptotic behavior of the confiuent hypergeometric function, 

F(a,7;^) « -f4-e^2"~'^, as l^l ^ 00, (61) 

the square integrability of the solutions (59) implies 

ivw_iv«_ r((ci-A„)/2)r(c2) 



ivf ivf r((c2-An)/2)r(ci)' 



(62) 



Now for the reference modes (9), we choose two eigenmodes belonging to an arbitrarily 
fixed energy E, given in terms of the solutions (59), as 



^(1) (^) ^ ^/ A 0(1) (1^1) [e(x) - e(-x)] , 

muj (g3) 



<p^'Kx) = ^—<l>^'\\x\). 



C2 - Ci 
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From F{a, 7; 2;) = 1 + 0{z) as 2; — > 0, the boundary vectors (10) are found to be 

which are finite despite that the solution (60) is divergent at the singularity. Moreover, 
we observe from the relations (62) and (64) that the two vectors and \1/ are again 
proportional to each other, and hence if we write = the boundary condition (2) 
becomes (32) as before. Combining (62) and (64), one finds 

^r((ci-An)/2)r(c2) ^gg^ 



C2-C1V n r((c2-A„)/2)r(ci)' 

and the spectrum {E^ — A^^cu} is determined by the same condition as in (37). Note 
that, again, the spectrum depends only on the two angle parameters (6'_|_, 9-). 

At some points of the angles, the spectrum becomes particularly simple. For instance, 
at {9+,9-) = (0,0) {i.e., U = /), we obtain En = {2n + 02)^100 and that the eigenstates are 
given by 0n (|ic|) both on the positive and negative half lines (and hence each level is doubly 
degenerate). Similarly, at (0+, 6'_) = (7r,7r) (i.e., U — —I), we find E^ = (2n + ci)hu) and 
that the eigenstates are (Z)!^"* (|a;|) which are also doubly degenerate. As mentioned earlier, 
this corresponds to the Friedrichs extension and has been conventionally considered for the 
quantization of the system since Calogero [19]. On the other hand, at (0_|_,0_) = (0,7r), 
then we have two series of eigenstates, one with N^"^ = Nj^^ = and the other with 
N^^ = N^^ = 0, with eigenvalues 

E^^^ = {2n + l + a)nw, E^^ = {2n + 1 - a)huj , n = 0,l,..., (66) 

respectively. In particular, in the limit ^ — > we have a —>■ 1/2, which shows that 
our system recovers the spectrum of a harmonic oscillator. A complete reduction to the 
harmonic oscillator system is realized by choosing U = ai, where the eigenstates become 
£-2/^/2 times the familiar Hermite polynomials (for a detailed discussion on the smooth 
limit to the harmonic oscillator, see [3]). In this respect, the extension provided by t/ = ui 
causes no obstacle at the singularity and is called 'the free case' in the analysis of point 
interactions. 
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4. Spectral space for parity invariant singular potentials 

The previous two examples share the property that the spectrum of the Hamiltonian 
is dependent only on the two parameters {9-^-, 9-) which are determined by the eigenvalues 
of the characteristic matrix U G U{2). This has been observed also for point interaction 
[18], and can be shown to hold for any singiilar potential V{x) characterized by U{2), as 
long as it is parity invariant V{—x) = V{x). Indeed, we have the following 

Theorem. If the Schrodinger operator H on X\{0} has a (measurable and locally inte- 
grable) parity invariant potential V{x) = V{—x), and is in the limit-circle case at x = 
and in the limit-point case for \x\ —>■ oo, then its spectrum on a self-adjoint domain Vjj is 
uniquely determined by the eigenvalues of the characteristic matrix U e C/(2). 

Proof. The proof is done simply by putting the argument of the examples in the general 
context. Let ipn £ -^^(^) be a normalizable solution of the Schrodinger equation (19) 
with eigenvalue which is subject to the boundary condition (2) specified by the matrix 
U. Let also {0n ^(a^)}s=i,2 be a fundamental system of real solutions for a; > with the 
same En- These solutions are not necessarily subject to (2), and are chosen to satisfy 
W[4)n\(l>n^] = 1- In terms of these, the general solution of V'n can be given in the 
form (60) because of the parity invariance of the potential, V{—x) = V{x). Since (at 
least one but generically both of) the basis solutions 0n ^ (x) become divergent as a; — > oo 
as dictated by the uniqueness of the solution at the limit-point infinity, one needs to 
arrange the coefficients, N^^ and so that the divergence of the two terms cancel 

each other in the limits x — > ±oo. From this one deduces the equality of the two ratios 
iv4^^ ivf = : ivi^\ that is, 

Arf = «ArW, Nl'^ = aNi'\ (67) 

with some a e M U {oo}. We also define our reference modes (9) as 

^W(x) = </>W(|x|)[e(x)-e(-x)], 
^(^)(x) = 0(^)(|x|), 

using two real eigenmodes 0*-^-* for a; > corresponding to an arbitrary eigenvalue E, 
satisfying = 1. Then, one obtains the following relations for the Wronskians: 
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With the help of these relations (69) and (67), one can compute the boundary vectors (10) 
to find 

One thus sees that the two complex boundary vectors are proportional to each other, 
= ^ ^f, with a constant ^ that is specified by (67) and (70). Once this is established, the 
rest of the argument is already given in the examples. Namely, the boundary condition (2) 
now becomes (32) and, hence, a nontrivial solution is obtained if (33) is fulfilled. In terms 
of the parameters (36) the spectrum condition reads (37). This proves the statement of the 
theorem for the bound states, since the diagonal slots of D are nothing but the eigenvalues 
of the matrix U . In the end, we recall the fact that the continuous spectrum is independent 
of t/, since all self-adjoint extensions of a symmetric operator admit the same continuous 
spectrum [7]. Q.E.D. 

We note that the space of possible spectra is therefore the space of the eigenvalues of 
[/, which is U{1) x [/(1)/Z2 (where Z2 is the factor of interchanging the two eigenvalues) 
forming a Mobius strip with boundary [18]. This theorem implies that, since the separated 
subfamily O3 contains all possible diagonal [/, the probability flow through the singularity 
plays no role as long as the spectra of parity invariant systems are concerned. In other 
words, the variety of the spectra is exhausted by systems consisting of two separated half 
lines, when all possible conditions at the boundary, i.e., the U{1) family of boundary 
conditions each, are allowed on both sides. 



16 



5. Summary and Discussions 



In this paper we presented the generahzed connection conditions for singular poten- 
tials characterized by the matrix U E U{2) in the form (2) with improved boundary vectors 
(10). An essential point in our generalization is the use of Wronskians in the boundary 
vectors (10) which are well-defined even in the limits to the singularity a; = 0, in contrast 
to the earlier ones (3) which become ill-defined in the limits. Using the generalized con- 
nection conditions, we examined two models, the one dimensional hydrogen atom and the 
harmonic oscillator with inverse square potential, which are solvable and yet so far have 
yielded conflicting results in the spectrum. Our analysis shows that the spectrum varies 
according to the connection conditions adopted, and that the possible spectra form a 2- 
parameter subspace (Mobius strip with boundary) in the entire U{2) family. The confusion 
on the spcctnmi is therefore resolved once we understand which connection conditions — 
if formulated in the form (2) — one is using in the analysis. We also note that in our 
connection conditions the parameters in the matrix U, when combined with the decom- 
position (34), bear direct physical meanings [18]. Indeed, we have already seen this in the 
spectral theorem in sect. 4 in that the two parameters in the diagonal piece D represent 
the two independent scales of the system. 

As a flnal remark, we wish to mention that the whole prescription for the connection 
conditions remains valid even for systems with a 'black box', not just a singular point. 
Namely, if there is a blank interval I = [— £0) ^o] with some small £o > on a line, then its 
quantum mechanical description can be given by means of our prescription if one replaces 
X\{0} with X\I, and the relevant Wronskians to be used in the boundary vectors are 
W^[0*)V']±eo instead of V'lio- This seemingly innocent modification has a practical 

consequence, since this allows us to introduce a rich U (2) structure to strong singularities 
that are in the limit-point case and are therefore originally essentially self-adjoint with no 
ambiguity in connection conditions. The appearance of the four-parameter U (2) freedom 
in choosing a singular potential may provide a useful theoretical framework for describing 
such singularities in quantum phenomena that seem to belong to the limit-point case by 
their behavior at intermediate length scales but have a richer structure at very short length 
scales. We mention that several recent approaches address the question of how to introduce 
nontrivial structure to limit-point singularities, see, e.g., [20]. 

Acknowledgement: I.T. is indebted to H. Miyazaki for useful comments. T.C. thanks 
members of the Theory Group of KEK for the hospitality extended to him during his stay. 
This work has been supported in part by the Grant-in-Aid for Scientific Research (C) 
(Nos. 10640301 and 13640413) and that on Priority Areas (No. 13135206) by the Japanese 
Ministry of Education, Science, Sports and Culture. 



17 



References 



[1] A. Fuhrer, S. Lusher, T. Ihn, T. Heinzel, K. Ensslin, W. Wegscheider and M. Bichler, 
Nature 413 (2001) 822. 

[2] T. Cheon, T. Fiilop and 1. Tsutsui, Ann. Phys. 294 (2001) 1. 

[3] H. Miyazaki and 1. Tsutsui, Ann. Phys. 299 (2002) 78. 

[4] A.N. Gordeyev and S.C. Chhajlany, J. Phys. A 30 (1997) 6893. 

[5] L. Lathouwers, J. Math. Phys. 16 (1975) 1393. 

[6] M. Reed and B. Simon, "Methods of Modern Mathematical Physics II, Fourier anal- 
ysis, self-adjointness" , Academic Press, New York, 1975. 

[7] N.I. Akhiezer and I.M. Glazman, "Theory of Linear Operators in Hilbert Space", 
Vol.11, Pitman Advanced Publishing Program, Boston, 1981. 

[8] G. Bonneau, J. Faraut and G. Valent, Self-adjoint extensions of operators and 
the teaching of quantum mechanics, LPTHE preprint PAR/LPTHE/99-43, |quant- 
I ph/0103T53 . 

[9] T. Fiilop and I. Tsutsui, Phys. Lett. A264 (2000) 366. 
[10] A.M. KraU, J. Differential Equations 45 (1982) 128. 

[11] S. Albeverio, F. Gesztesy, R. H0egh-Krohn and H. Holden, "Solvable Models in Quan- 
tum Mechanics" , Springer, New York, 1988. 

[12] F. Relhch, Math. Ann. 122 (1951) 343. 

[13] C.T. Fulton, Trans. Amer. Math. Soc. 229 (1977) 51; Proc. Roy. Soc. Eddinghurgh 
Sect A 87 (1980) 1. 

[14] A.N. Kochubei, Siberian Math. J. 32 (1991) 401. 

[15] R. Loudon, Amer. J. Phys. 27 (1959) 649. 

[16] W. Fisher, H. Leschke and P. Miiller, J. Math. Phys. 36 (1995) 2313. 

[17] F. Gesztesy, J. Phys. A 13 (1980) 867. 

[18] I. Tsutsui, T. Fiilop and T. Cheon, J. Math. Phys. 42 (2001) 5687. 

[19] F. Calogero, J. Math. Phys. 10 (1969) 2191, 2197; 12 (1971) 419. 

[20] P. Kurasov and K. Watanabe, On Rank One H -^-Perturbations of Positive Self- 
adjoint Operators, Stockholm Univ. Res. Rep. Math. 1999/10; On H-^-perturbations 
of self- adjoint operators, Stockholm Univ. Res. Rep. Math. 2000/16. 



18 



